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We present a study of the lensing properties of two-dimensional (2-D) photonic 

quasicrystal (PQC) slabs made of dielectric cylinders arranged according to a 12- 

O fold-symmetric square-triangle aperiodic tiling. Our full-wave numerical analysis 

^ confirms the results recently emerged in the technical literature and, in particular, 

the possibility of achieving focusing effects within several frequency regions. How- 

pi I ever, contrary to the original interpretation, such focusing effects turn out to be 

critically associated to local symmetry points in the PQC slab, and strongly depen- 

psl dent on its thickness and termination. Nevertheless, our study reveals the presence 

^ of some peculiar properties, like the ability to focus the light even for slabs with 

~f^ a reduced lateral width, or beaming effects, which render PQC slabs potentially 

Cn| interesting and worth of deeper investigation. 
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a 1 Introduction 



Photonic crystals (PCs), featuring a periodically modulated refractive index, 
represent an active research area since nearly two decades [Tp] . The strong 
interest toward these structures is motivated by their capability of inhibiting 



* This paper was presented at the PECS VII Conference. Corresponding author: 
E. Di Gennaro, Tel. -^39-081-7682661, Fax -^39-081-2391821. 
Email address: emiliano@na.infn.it (E. Di Gennaro). 

Preprint suhmitted to Elsevier 5 February 2008 



the spontaneous emission of light, which potentially renders them the building- 
block components of an all-optical circuitry [3] . 

In PCs, multiple wave scattering at frequencies near the Bragg condition pre- 
vents propagation in certain directions, producing a bandgap. Anomalous (even 
negative) refraction and subwavelength focusing ( "superlensing" ) properties 
for PCs have been theoretically predicted (see, e.g., P|l5|l6|f7] ) and experimen- 
tally verified (cf. [SE])- The theoretical foundations of the (analytic and nu- 
merical) modeling, prediction, and physical understanding of the above phe- 
nomena are based on well-established tools and concepts developed for the 
study of wave dynamics in periodic structures, such as Bloch theorem, unit 
cell, Brillouin zone, equifrequency surfaces (EPS), etc. 

With specific reference to lensing applications, two different approaches have 
been presented to obtain subwavelength resolution using a PC slab. The first 
one assumes that, under suitable conditions (single beam propagation, nega- 
tive group velocity in the specific frequency band, and circular EPS), a PC can 
behave like a homogeneous material with a negative refractive index n = —1. 
These features are likely found near a frequency band edge and in the case 
of high dielectric contrast [1]. Purthermore, modifying the slab terminations, 
the impedance mismatch with the surrounding medium can be minimized 
[To] and/or surface waves can be excited allowing (partial) reconstruction of 
evanescent waves [11]. Such a slab is not a perfect lens [12], in the sense that 
the image has a finite resolution, since not all evanescent components can be 
restored. Nevertheless, the focus position follows the simple ray-optical con- 
struction as for a fiat lens with n = —1 [GflS] . The second approach relaxes 
some of the previous requirements (negative group velocity and isotropic re- 
fraction properties) and achieves "all-angle negative refraction" (AANR) with- 
out an effective negative index, provided that the EPS are all convex and larger 
than the frequency contour pertaining to the surrounding medium [7]. This 
condition is still sufficient for a PC slab to focus with subwavelength resolu- 
tion, if the termination is designed so as to sustain surface modes. It is worth 
noticing that, in this case, the focus position does not follow the ray-optical 
construction [T3] and is restricted 



Recently, there has been an increasing interest toward the study of aperi- 
odically ordered PCs. Such structures are inspired by the "quasicrystals" in 
solid-state physics [T5|IT6] and the related theory of "aperiodic tilings" [T7] . 
and are accordingly referred to as "photonic quasicrystals" (PQCs). Like their 
solid-state physics counterparts, PQC structures can exhibit long-range (e.g., 
quasiperiodic) order and weak (e.g., statistical and/or local) rotational sym- 
metries which are not bounded by the crystallographic restriction and can 
thus be considerably higher than those achievable in periodic PCs p[8]. Ac- 
cordingly, it has been proposed to utilize PQC structures in order to maintain 
the periodic-like scattering of light, while introducing additional geometric 



degrees of freedom potentially useful for performance control and optimiza- 
tion. However, it must be said that the lack of periodicity renders the study 
of PQCs very complex and computationally-demanding. Although some con- 
cepts developed for periodic PCs can be used for the analysis of the properties 
of PQCs [T9|20] . a rigorous extension of a Bloch-type theorem (and associated 
tools and concepts) does not exist. 

Recent studies of the EM properties of 8-fold ^21], 10-fold [T91I22] and 12-fold 
[IH] PQCs have shown that a gap in the density of states, accompanied by sig- 
nificant dips in the transmission spectra, can be achieved. PQCs also allow, in 
principle, for a higher degree of flexibility and tunability for defect-type [23)121] 
and intrinsic [25] localized modes. Moreover, examples of subwavelength fo- 
cusing properties have been presented for 2-D PQCs exhibiting various (8-fold, 
10- fold, and 12-fold) rotational symmetries [^UfTT] . Similar results were also 
obtained in the case of acoustic waves 



Although the above studies on the focusing properties of PQC slabs reveal 
some appealing properties, such as polarization-insensitive and non-near-field 
imaging, we believe that the original phenomenology interpretation provided 
in [2n], within the framework of "effective negative refractive- index and evanes- 
cent wave amplification," is questionable and deserves a deeper investigation. 
In this paper, we present some representative results from a comprehensive 
parametric study of the focusing properties of a 12-fold symmetric PQC slabs 
(as in [26]), based on full- wave numerical simulations (Fourier-Bessel mul- 
tipolar expansion). In particular, with reference to the line-source imaging 
scenario, we investigate the critical role of the global vs. local symmetry of 
the lattice (which is merely glossed over in [27]), by varying the PQC slab 
thickness, lateral width, and termination, as well as the source position in 
both parallel and orthogonal directions with respect to the slab surface. 



2 Samples and Methods 



As in [26], the 2-D PQC of interest is made of infinitely- long dielectric cylinders 
with relative permittivity £r = 8.6 and radius r = 0.3a placed (in vacuum) 
at the vertices of a 12-fold-symmetric aperiodic tiling with lattice constant 
(tile side-length) a, generated by the so-called Stampfli inflation rules [29] 
illustrated in Fig. Ilia). The (recursive) generation algorithm is based on a 
parent tiling represented by the central gray-shaded dodecagon. Scaling up 
this structure by an inflation factor u = -\/3 + 2, one obtains a big parent (red 
dashed dodecagon in Fig. [l|a)). Copies of the original tiling are subsequently 
placed at each vertex of the big parent, and the process is iterated up to the 
desired tiling extension. From a suitably-sized tiling, we extract several PQC 
slabs of different sizes, with a thickness t (along the x-axis) ranging between 





Fig. 1. Dodecagonal PQC geometry, (a): Illustration of the Stampfli inflation rule. 
The gray-shaded dodecagon in the center represents the parent tiling from which a 
big parent (red dashed lines) is generated by applying a scaling factor v = \/3 + 2. 
(b): Two examples of PQC slabs of different thickness (black-dotted and red-dashed 
rectangles) extracted from the tiling. The green- full- dots dodecagon in the center 
corresponds to the parent tiling in (a). 

~ 7a and ~ 11a (see the two rectangles in Fig. [lib), corresponding to the two 
examples presented in [26]) and a lateral width h (along the y-axis) varying 
between ~ 11a and ~ 30a. Our numerical simulations are carried out via a 2-D 
full-wave method based on a multipolar expansion of the fields around (and 
inside) each cylinder [5U] . This method is closely related to the Korringa-Kohn- 
Rostocker method used in solid-state physics |3l] (which is itself derived from 
a pioneering work by Lord Rayleigh in electrostatics p2]), and has also been 
recently applied to the study of the local density of states in finite-size PCs 
[33,34j and PQCs [22] • In the present investigation, only E-type polarization 
is considered, i.e. the electric field is assumed to be parallel to the cylinders. 



3 Representative Results 



As a preliminary step, we studied the focusing properties of the two PQC 
slabs in Fig. ^[b), for a fixed line-source position, and varying the normalized 
frequency a/A (with A being the vacuum wavelength) within a range between 
the first and second bandgaps (0.354 < a/A < 0.533). We found several fre- 
quency regions (including those in ^26j) where a clear focus was visible. We 
note that this feature represents a first remarkable difference with the PC case, 
where the focusing regions tend to be more rare and well separated. In PQCs, 
likely due to their inherent self-similar nature [20], focusing effects seem to 
occur more densely in the spectrum. In all simulations below, attention is fo- 
cused on the normalized frequency a/A = 0.394, corresponding to the value for 
which a refractive behavior emulating a homogeneous medium with n = — 1 
was estimated in 
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Fig. 2. Intensity field maps for various symmetric slab configurations with different 
values of the thickness t. (a): t = 11a; (b): t = 9.8a; (c): t = 9a; (d): t = 8.8a; 
(e): t = 8a; (f): t = 7a. An electric line-source is placed at y = and a distance 
dx = t/2 from the vacuum-slab interface. The crosses mark the source position and 
the ray-optical (n = —1) prediction for the focus location. The color-scale limits 
in the plots are chosen so as to properly reveal details in the internal and external 
fields that would be otherwise overwhelmed by the dominant features. 

3.1 Effects of slab thickness 



In order to gain some insight into the focusing properties, we then consid- 
ered PQC slabs of different thickness. The line-source was always placed at a 
distance dx = t/2 from the vacuum-slab interface, and in correspondence of 
the center of the parent tiling [y = 0). From a simple ray-optical construc- 
tion (for fiat lenses with n = —1), the focus is accordingly expected to be 
located in a symmetric position with respect to the slab. The slab thickness 
was changed by removing cylinders in such a way to preserve the inversion 
symmetry with respect to the y-axis, thereby obtaining six different config- 
urations, which include as extremes the two samples shown in Fig. fltb). As 
a first result, our simulations confirmed the results in [2S], as shown in the 
(intensity) field maps in Figs, ^a) and|2](f). A focus is clearly observable at a 
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Fig. 3. As in Fig. [2] but for asymmetric configurations, (a): t 
(c): t = 9a; (d): t = 8a; (e): t = 7a; (f): t = 6.6a. 



10.3a; (b): t = 9.8a; 



position very close to the ray-optical prediction (marked by a cross symbol). 
However, for intermediate thickness (see Figs. 

[2Fb)-[2](e)), the focusing properties become rather questionable. Moreover, in 
all the simulations, the reflection coefficient in the source plane is considerably 
nonzero, indicating the presence of a significant impedance mismatch at the 
vacuum-slab interface. Furthermore, the field distributions inside the slab turn 
out to be significantly different, suggesting a strong dependence from the whole 
slab configuration, including its termination. In this connection, we note that it 
is not possible to change the PQC lens thickness without changing at least one 
termination. In Figs.|2](c) and[2|d), a feature that resembles a "beaming effect" 
can be observed, together with a strong field intensity along the vacuum-slab 
interface on the image side. In periodic PC slabs, such a phenomenon has 
been related to the presence of localized photon states on the slab surface [7] , 
stemming from an overall resonance in the structure. 

In Fig. |3} the field maps pertaining to other six different PQC slabs of various 
thickness are displayed. In these examples, the slab thickness was changed 
by removing dielectric cylinders only from the image side, thereby destroying 



the inversion symmetry with respect to the y-axis. Also in this case, the field 
distribution on the image side appears to depend critically on the structure 
geometry. Although a clear single focus cannot be identified, several spots 
(Figs, jsj^a) and[3|^d)), or beaming effects (Figs. |3](b), [3|e), andjsj^f)) are visi- 
ble, indicating that the inversion symmetry does play a key role in the focus 
formation. 



3. 2 Effects of source position 



We then fixed the slab geometry to that of Fig. Wia) (thicker slab), that was 
found to exhibit focusing properties, and changed the source position. In case 
of a fiat lens with n = —1, simple ray-optics predicts the source- image dis- 
tance to remain constant and equal to twice the lens thickness [6] . Results for 
periodic PCs showing an almost isotropic refractive response (see, e.g., pj]) 
agree fairly well with this prediction. First, we kept the source in a symmetric 
position along the transverse direction {y = 0), varying the distance dx from 
the slab from 0.75a up to 10.5a. In Fig. |4], some representative field maps 
are shown. A variety of effects is observed, ranging from the absence of focus 
(Fig. |4]^a)) to multiple spots (Figs. |4]^c) and|4](d)), but always in the near- 
field region, raising further concerns about the interpretation in [26]. Similar 
results [13] were also observed for PCs, in the case of an anisotropic regime 
with positive or negative refraction properties, where the near-field focusing 
was found to be dominated by a self-coUimation effect. 

Next, for a fixed source-slab distance d^ = t/2, we moved the source parallel 
to the vacuum-slab interface. Representative results are shown in Fig. [5j In 
PCs, for source displacements that preserve the distance from the slab and the 
structure periodicity, the focus remains unaffected, and its position follows the 
source location. Conversely, in our PQC case, when the x-axis symmetry is 
broken, the focus undergoes a rapid deterioration, and is completely destroyed 
for a source displacement of dy = Qa (see Fig. |5](d)). These findings highlight 
the importance of keeping the source on the same axis of the central symmetry 
point in the 12-fold tiling (highlighted in Fig.[l|b)), as also observed in [27]. For 
the thicker slab configuration in Fig. W(h), we also performed some simulations 
involving a coUimated Gaussian beam impinging on three different points of 
the PQC slab interface with three different incidence angles. Results shown in 
Fig. [6} displaying complex multi-beam features in the transmitted field that are 
strongly dependent on the incidence point, highlight the absence of a clear-cut 
refractive behavior, and explain the strong dependence of the focusing effects 
on the source distance from the slab (cf. Fig. Ill) and lateral displacement (cf. 
Fig.|5|. 
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Fig. 4. As in Fig. ^a) (i.e., thicker slab, cf. red dashed rectangle in Fig. U[b)), but 
for various source distances dx from the vacuum-slab interface, (a): dx = 8.3a; (b): 
dx = 7.1a; (c): dx = 3.6a; (d): dx = 2.5a. The source y-position is maintained at 
y = 0. 
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Fig. 5. As in Fig. El but for various source displacements dy along the y-axis. (a) 
dy = 0.75a; (b): dy = 1.5a; (c): dy = 3a; (d): dy = 6a. The source distance from the 
vacuum-slab interface is maintained at dx = 5.5a (i.e., half the slab thickness). 
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Fig. 6. As in Fig.^a), but for Gaussian-beam incidence at various lateral positions 
m at the slab interface, and various incidence angles 6i (with respect to the rr-axis). 
The Gaussian beam minimum spot-size is ~ 2.2a and the incidence point is located 
at the Rayleigh distance, (a): yi = 0, 0, = 10°; (b): yi = -l.Ua, Oi = 10°; (c) 
Vi = -4.88a, Oi = 10°; (d): y, = 0, 0^ = 45°; (e): yi = -2.44a, e, = 45°; (f) 
yi = -4.88a, Oi = 45°; (g): y, = 0, 0^ = 60°; (h): y, = -2.44a, Oi = 60°; (i) 
y. = -4.88a, Oi = 60°. 



3. 3 Effects of slab lateral width 



To better illustrate the role of the local symmetry point, we also varied the slab 
lateral width asymmetrically (i.e., removing cylinders only from one side of the 
l/-axis) while keeping the source in the central position [d^ = t/2 and y = 0). 
Results are shown in Fig. [7j In this case, the focus turns out to be rather 
stable, in spite of the x-symmetry breaking. These results clearly reveal the 
key role played in the focus formation by short-range interactions involving 
a neighborhood of the parent tiling (highlighted via green full dots in Fig. 
fiyb)). This can be also observed considering a configuration where the slab 
lateral width is symmetrically reduced to ~ 11a with respect to the y axis, as 
shown in Fig. [8J Again, in spite of the considerably reduced lateral width, the 
focusing properties associated to the symmetry local point located at (x = 
0,y = 0) appear to be rather robust, as if their occurrence was restricted 
to a hmited range of incidence angles. This renders these PQC structures 
particularly appealing for the development of compact optical systems. 




Fig. 7. As in Fig. ^a), but for various asymmetric slab configurations with different 
lateral width h. (a): h = 29.4a; (b): h = 26.7a; (c): h = 23.9a; (d): h = 20.6a. 




Fig. 8. As in Fig. [7j but reducing the lateral width to ~ 11a symmetrically with 
respect to the y-axis. 

4 Conclusions 



To sum up, a comprehensive parametric numerical study of the lensing prop- 
erties of quasiperiodic dodecagonal PQC slabs has been carried out. Results 
confirm that these structures possess focusing properties in several frequency 
regions, even if critically associated to local symmetry points (and neigh- 
boring regions) existing in the tiling lattice. Indeed, our full-wave numerical 
simulations show that varying the slab thickness or the source position, the 
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field distributions inside the slab and in the image plane change dramatically. 
The original interpretation in [2S] , within the framework of "effective negative 
refractive-index and evanescent wave amplification," does not appear adequate 
to describe the underlying phenomenologies, which entail complex near-field 
scattering effects and short-range interactions. Nevertheless, this particular 
PQC shows some peculiarities in the EM response that deserve further inves- 
tigation. In particular, applications to lenses with small aperture or antennas 
with highly directional beaming could be envisaged. 

Also of interest it is the parametric study of the refraction and focusing prop- 
erties of other classes of PQCs based on periodic (e.g., Archimedean [35]) and 
aperiodic (e.g., Penrose p^|l22] . octagonal [2T], circular [36]) tilings character- 
ized by diverse degrees of (local and/or statistical) rotational symmetry. 
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